A coupled order parameter system on a scale-free network 
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The system of two scalar order parameters on a complex scale-free network is analyzed in the 
spirit of Landau theory. To add a microscopic background to the phenomenological approach we 
also study a particular spin Hamiltonian that leads to coupled scalar order behavior using the mean 
field approximation. Our results show that the system is characterized by either of two types of 
ordering: either one of the two order parameters is zero or both are non-zero but have the same 
value. While the critical exponents do not differ from those of a model with a single order parameter 
on a scale free network there are notable differences for the amplitude ratios and susceptibilities. 
Another peculiarity of the model is that the transverse susceptibility is divergent at all T < T c , 
when 0(n) symmetry is present. This behavior is related to the appearance of Goldstone modes. 
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I. INTRODUCTION 

The topology of many natural and man-made networks 
(social networks, biological, technological and trans- 
portation systems) strongly differs from the topology of 
regular lattices or even random graphs. Often these net- 
works show scale-free behavior [l], 0, Q - the probability 
of a randomly chosen node to have a degree fc (to have fc 
links) follows a power law 



P(k) = Ak 



-A 



(1) 



Other integral parts of many real networks are a small- 
world effect and high clustering, resulting in specific fea- 
tures of cooperative phenomena on such systems. This 
has sparked interest in the analysis of different spin mod- 
els on complex networks [4|. Moreover, such models have 
interesting applications. For example, the opinion of each 
individual of a social network may be represented by an 
Ising spin e.g. for simple YES or NO alternatives. Such 
a model may describe phenomenon of opinion formation, 
namely, individuals change their mind under the influ- 
ence of their acquaintances [5|. 

The properties of the order-disorder phase transition 
of the Ising model on complex networks strongly depend 
on the node degree distribution ([1]). Numerical simula- 
tions @ and analytical calculations Q of the Ising model 
on Barabasi- Albert scale-free networks (A = 3) as well as 
different analytical approaches [1, [9[ and Monte Carlo 
simulations [1 01 ] for the Ising model on networks with ar- 
bitrary degree distributions have been performed. Three 
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types of behavior were found depending on the respec- 
tive behavior of the moments (k 2 ) and (fc 4 ) of the degree 
distribution, which are related to the value of the A ex- 
ponent. Namely, if (fc 4 ) and (fc 2 ) are finite (A > 5), 
the behavior of the system is described by the standard 
mean-field critical exponents. If (fc 4 ) diverges and (fc 2 ) is 
finite, which corresponds to 3 < A < 5, the critical behav- 
ior is governed by either mean fields exponents with loga- 
rithmic corrections (A = 5) or by nontrivial A-dependent 
critical exponents. Finally, if both (fc 4 ) and (fc 2 ) diverge 
(2 < A < 3), the critical temperature becomes divergent 
(for the infinite-size networks) and the system is always 
ordered. Furthermore, other models on scale free net- 
works, namely the XY [ll[ and Potts [HI, 0j| models (for 
a more detailed list see e.g. also show peculiarities 

depending on the value of A. 

Rather recently, critical phenomena on complex net- 
works have been studied in the spirit of Landau theory 
[l5j . The power of the latter is that it is independent of 
the origin of the interactions between the particles, and 
therefore it may be applied to a wide range of systems. 
The main new feature of the phenomenological theory 
of critical phenomena on complex networks, that differs 
from the standard Landau theory is the dependence of 
the coefficients on the moments (fc 1 ) of the degree distri- 
bution IJ). 

Landau theory for two interacting scalar order param- 
eters is widely used to analyze systems with several pos- 
sible types of ordering (e.g. ferromagnetic and antiferro- 
magnetic, ferroelectric and ferromagnetic, or structural 
and magnetic ordering). Such combination of order pa- 
rameters may be described by a model of two scalar order 
parameters x\,xi, which are coupled Assuming 
that the Landau free energy is analytic and symmetric 
with respect to the signs of x\ and X2 the lowest order 
coupling is biquadratic 

$(f,T) = %{T - T c )\x\ 2 + h -\xf + -x\x\. (2) 
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where x = (x%,X2), \x\ 2 = x\ + x?,, a, b, c are phenomeno- 
logical Landau parameters, T and T c are temperature 
and critical temperature, correspondingly. A possible 
application of this model of coupled order parameters 
on a social network may reflect the coupling between the 
preferences for a candidate and a party in an election (or 
similar scenarios) [Hj]. The free energy @ corresponds 
to the free energy of an n- vector anisotropic cubic model 
in the case n = 2. The latter is obtained from the 0(n) 
invariant free energy by adding invariants of the symme- 
try group B n of the n-dimensional hypercube [l9j . 

The aim of our work is to generalize the Landau theory 
for models on complex networks [l5j to the case of two 
interacting order parameters with a free energy symme- 
try given by ([2]). The structure of our paper is as follows. 
The next section ((11]) lays out the basic assumptions of 
the theory and the peculiarities of the free energy con- 
struction and compares the approach with a microscopic 
model. Section IIIII describes the stable states and the 
phase diagrams of the system. The behavior of the ther- 
modynamic functions, isothermal susceptibilities and the 
heat capacity is described in section IIVI We conclude 
with an outlook in section [V] Some details of our calcu- 
lations are given in appendices and [B] 



II. FREE ENERGY 

This section is devoted to the construction of a gen- 
eralized Landau theory for a system with two coupled 
order parameters on a network pi A[) . Besides, we derive 
a corresponding free energy starting with a microscopic 
spin Hamiltonian and compare both approaches pi Bp . 



A. Generalized Landau theory 

In the spirit of the Landau approach we assume that 
the system may display some ordering which can be 
quantitatively characterized by two order parameters 
x\ and X2- For convenience let us introduce a vector 
x = (xx,X2)- Following the work of Ref. we assume 
that the Landau free energy is not only a function of the 
order parameters, conjugated field h, and temperature 
but also depends on the node degree distribution P(k) 



$(x,T,h) 



dkP(k)f(x,kx) - hx, (3) 



where f(x, kx) represents the contribution to the free 
energy of an individual node of degree k, and k max is 
the maximal node degree of the network. Note, that 
k m ax — > oo is implied for an infinite size system with a 
power-law node degree distribution as in Eq. (Q]). That 
f(x, kx) depends not only on the order parameters X\ 
and X2, but also on kx may be understood by simple 
reasoning. It reflects that any node with k neighbors is 
subject to a field kx of these neighbors. 



The next basic assumption in the case of a scalar order 
parameter x is that f(x, kx) is an analytic function of x 
and kx [HI]. In the case of two order parameters we 
assume that f(x,kx) is now an analytic function of Xi, 
X2, kx\ and kx2 and may be represented as a series in 
their powers 



f(x , kx) 



flil2m 1 m 2 x l x 2 2 1 {kX2) 



li : l2,mi : nt2—0 



(4) 

where fi 1 i 2 m 1 m 2 are functions which in general may de- 
pend on the temperature T and an external field h. More- 
over, some relations between these coefficients are im- 
plied by the symmetry of the system as described by 
equation In this case the function f(x,kx) may be 
represented as 



f(x,kx) = f Q + ^2a,k l \x\ 2 + ^2b. t k 



l, x\ 4 



(5) 



i=0 
4 2 



i=0 



where /o, a^, bi, Ci are convenient notations for the coef- 
ficients 

Naturally, the free energy ((3]) must be finite if the or- 
der parameters are finite. This condition is satisfied in 
particular if the behavior of the function f(x, kx) at large 
k\x\ — > oo is bounded by 



/(a?, kx) ~ k\x\, k\x\ 



oo. 



(6) 



The assumptions (J3]) and © as well as condition ^ 
serve as our basis to analyze the phase transitions in the 
coupled order parameter system following the standard 
approach of Landau theory 20] . 

Substituting ([5]) into ([3]) and taking into account that 
the coefficient of \x\ 2 in the free energy is equal to zero 
at the critical point, the equation for the critical tem- 
perature T c as function of the moments of the degree 
distribution is found in the same manner as in the case 
of a scalar order parameter [l5| to be 



a (T c ) + ai(T c )(fc) + a 2 (T c )(k 2 ) = 0. 



(7) 



If ao(T c ) — 0, the critical temperature is a function of 
(k 2 )/(k). This statement is in accordance with the ex- 
act result for the Ising model on networks obtained ana- 
lytically n, n and confirmed numerically [10j , where T c 
follows 



1 



111 



(k 2 



(k 2 ) - 2(k) 



(8) 



Before we embark to calculate the free energy let us 
discuss an essential point that is the origin of many of 
the peculiarities of cooperative phenomena on networks. 
For scale- free networks with a node degree distribution as 
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in Eq. (p} one finds in general that all moments (k l ) with 
i < X — 1 are finite, whereas all moments with i > X — 1 
diverge. If we restrict the series in Eq. (f5]) to the fourth 
power of the order parameter, there are no relevant di- 
vergent moments for A > 5. Nevertheless, if (k A ) or lower 
moments of the degree distribution are divergent (A < 5), 
as often found for real networks, the free energy © at 
the first sight may seem to be infinite for any nonzero val- 
ues of the order parameters, a behavior which certainly is 
unphysical. In fact, the correct way to calculate the free 
energy is to take into account all the orders of the series 
([5]). This procedure ensures to a behavior of the func- 
tion f{x, kx) at large values of k\x\ — > oo as described by 
equation ((6]). Therefore, we collect all terms in equation 
([5|) containing k % with i > A — 1 together with the highest 
orders of the series (JS|) in a function g(x, kx): 



12 



f(x,kx) = fo + y ai¥\xY + 



i=0 



14 

E 

i=0 



hk l 



(9) 



+ E CifcJ E x M + 9(x,kx). 

Here i 2 is the maximal integer that satisfies both con- 
ditions i 2 < 2 and i 2 < A — 1. Respectively, i 4 is the 
maximal integer that satisfies both i 4 < 4 and i 4 < A — 1. 
Now it is straight forward easy to integrate the part of 
f{x,kx) that does not include g(x,kx). Any peculiari- 
ties are connected with the integration of g(x,kx). Let 
us therefore investigate the properties of this function. 
Comparing © with §5§ one finds that for small values of 
k\x\ and for 3 < A < 5 this function behaves as 

g (x,kx) = b4k\x\) A +c 4 ^^-(k\x\) 4 , fc|*l->o. (10) 

For 2 < A < 3 one finds the following behavior: 

g(x,kx) = a 2 (k\x\) 2 + (b 2 + c 2 ^±^)\x\ 2 (k\x\) 2 



(63 



c 3 - 



?|(fc|x|) 3 



(11) 



(fT0| - (TT2|) ). Let us pass to a new variable y = k\x\, which 
ranges from |x| to infinity for infinite size networks. For 
a network with a power law node degree distribution (|TJ) 
one may then write 



dkP(k)g(x,k\x\) = A\x 



A-l 



^g{x,y). (13) 



As only the asymptotics of g{x, y) are fixed, let us write 

|S| dy ,„ s , s 

-jg&v), (14) 



\s\ v 



%9{x, y) 



%a{x, y) 



y 



y 



where e is a small positive number < e < \x\; obviously, 
both sides of this expression do not depend on e. The 
first term on its right-hand side is convergent, due to the 
asymptotic behavior (|12p . In the second term g(x, y) may 
be replaced by its expansion for small values of y (110p . 

CD- 

Following this procedure, one may obtain the free en- 
ergy. Details of the integration of the expression ([T4]) 
for the case 4 < A < 5 are given in Appendix [5] (ana- 
logue calculations can be performed for other values of 
A). In the following, we will consider zero external mag- 
netic field h = and drop the explicit ^.-dependence from 
our notations. Let us present the resulting expressions 
for the Landau free energy for different ranges of val- 
ues of A. We treat the cases (a) A > 5, (b) X = 5, (c) 
3 < A < 5, (d) X = 3 and (e) 2 < A < 3. As we will 
see, differences between usual Landau theory and that 
on a scale-free network become apparent starting from 
the marginal case A = 5. 

(a) Case A > 5 

In this case the free energy may be found easily by 
substituting ([5]) into ([3]) and performing the integration. 
The free energy reads: 



$(f,T) = f + -(T-T c )\x\' + — |x| 4 + — x{x l 2 . (15) 

The specific network properties are expressed by the co- 
efficients: 



-(T-T c ) = ai (k)+a 2 (k 2 ) 



(16) 



(64 



c 4 - 



%|£|) 4 , fc|f|-0. 



We do not consider the case A < 2 here as far as (k) is 
not defined. 

In order to satisfy condition © for a finite free en- 
ergy, the behavior of g(x, kx) is restricted for large val- 
ues of k\x\ — > 00 by the highest explicitly written term 
of f(x,kx) in Namely, for k\x\ — > 00 g{x,kx) is 

restricted by 



g{x, kx) 



(fc|f|) 3 , 4 < A < 5 
(k\x 



3 < A < 4 



(k\x\), 2<A<3. 



(12) 



To perform the integration of g{x, kx) in ((3]), note that it 
actually depends on x and k\x\, g(x, kx) = g(x,k\x\) (see 



b ix 1 = 4b 4 (k 4 ), c (x) =4c 4 (k 4 ). 



(17) 



As seen below, (fTHj) also holds for 3 < A < 5. 

(6) Case A = 5 

In this case the free energy reads: 

n M A ) 1 

$(f,T) = /o + -(T~T c )|f| 2 + — |f | 4 In - 



4 1 1 



x\x% In 



1 

W\ 



(18) 



In this marginal case the free energy displays logarithmic 
corrections to the standard mean-field behavior. The co- 
efficient a(T — T c ) is described by p^]) and the other 
coefficients are as follows 



b (X ) = 4,4(64 + C4), c 



(A) 



(19) 
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(c) Case 3 < A < 5 

Here, the free energy reads: 



6 (A) 



•I A — 1 



I A — 1 



$(x,T) = f +~(T-T c 

(20) 

In this case the free energy (|2T)]) explicitly depends on 
A. The coefficient a(T — T c ) is also described by (fl"6|) . 
whereas to get expressions for and from the in- 
tegration of g(x,kx), one needs to perform explicit cal- 
culations in parallel to those, presented in the Appendix 

m 



(cf) Case A = 3 

Here, the free energy reads: 

$(x,T) = f + C\x\ 2 -D\ocf 



E- 



(21) 



(e) Case 2 < A < 3 

In this case we find a free energy of the form 



$(f,T) = f Q + C 



D If 



A-l 



■ E 



}^\x\ x - 



(22) 



For the cases (d) and (e) we give explicitly only the 
expressions for D (for A = 3) and C (for 2 < A < 3): 



D 



-Aai 



C =a x {k) 



Aa 2 
3- A' 



(23) 



For an example how to calculate the other coefficients 
from the integration of g(x, kx), see Appendix 

Note, that for 2 < A < 3, the term of order \x\ 2 is no 
more the leading one. Terms of lower order of magnitude 
become relevant. In particular there is a term \x\ 2 In |a?| _1 
for A = 3, and a term |x| A_1 for 2 < A < 3. 

Before passing to the details of the phase diagram that 
results from the expressions for the Landau free energy 
([15]). (|i"8|). 1(20"] ) . l[2i" |) . (|22|) . we first proceed to show that 
the Landau free energy may also be derived from a spin 
system on a network by calculating its partition function 
in the simplest of approximations. 



B. Anisotropic Hamiltonian 

One of the ways to get the Landau free energy with 
two coupled scalar order parameters is to start from two 
coupled spin subsystems [211] . Another way is to consider 
a single spin system with a cubic anisotropy term. Let 
us here use the second option, considering a spin model 
on a complex network described by a Hamiltonian with 
an anisotropic term 



N 



(id) 



EE' 



(24) 



the notation Ylu j) denotes the summation over all pairs 
of connected nodes, the index v numbers the components 
of the two-component vector, si ■ Sj — J2v=i s f,i s v,j is 
a scalar product. Again, as above we will consider the 
case when the network node degree distribution obeys 
a power law decay (Q]). Note, that the Hamiltonian (j2"4")) 
represents an n- vector anisotropic cubic model [l9j in the 
case n = 2. 

Here, we consider the Hamiltonian (f2~4"|) in the spirit of 
a mean-field approach. Applying the mean field approach 
to a model that is defined on a regular lattice (equal 
degree k for all nodes), each node is characterized by the 
same mean spin (s) and experiences the effective field 
(k)(s) of its (k) neighbors. In the case of a complex 
network, this assumption may be applied only to nodes 
with the same degree: in the simplest approximation each 
fc-degree node experiences the same mean spin (s)k- In 
turn, the mean spin value per node (s) may be expressed 
in terms of (s)k as 



( S -)=£P(fc)< S -) 



k ■ 



(25) 



On the other hand, it can be found from the thermody- 
namical definition 



( dG(T,hy 
V dh >t 



(26) 



Here G(T, h) is the appropriate thermodynamical poten- 
tial and h is an external field. 

Node i experiences the effective field of its ki neigh- 
bors. This field may be quantitatively characterized by 
the mean value a % of the spins surrounding the i-th node 



(27) 



(J) 



Here the sum over j spans the ki nearest neighbors of 
node i. Now, in the spirit of the mean field theory one 
assumes that a 1 does not depend on the node number i 



1...N. 



(28) 



Note, that the above defined value a differs from the 
mean spin value per node (s). Equation (|26[) gives us the 
relation between the mean spin (s) and the effective spin 
a per neighbor and may be treated as a self-consistency 
equation. 

To proceed with the Hamiltonian ([24]) . we introduce 
the deviation of every spin component s v ^ from the cor- 
responding component of average spin per neighbor <?: 



As,,, = 



(29) 



where si and s"} are spins on nodes i and j correspond- Substituting s v i — er„ + As^j into the scalar product 
ingly, J and u are the coupling and anisotropy constants, in ([M)) and neglecting the terms of order (3((As) 2 ) we 
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arrive at the mean-field Hamiltonian: 

N 



i = l 



with 
H 



AIF 



j 2 2 

-J(k)a 2 - Jki^^UySv^ + u f 



(30) 



(31) 



Here, ct 2 = ct 2 + a 2 . Now, the partition function is re- 



duced to a product of single-site traces: 

N 

Zmf = II ^ e' H "^ T . 



(32) 



Here, the single-site trace Tr^(. . .) denotes the integration 
over all possible directions of si: 



last terms in (|36|0 - a property postulated in the Landau 
approach (see Section Hi A| . 

We now replace the sum over nodes in |35|) by a sum 
over node degrees 

N k max 

F & T ^nT, M k ^ = E PWffr ka). (39) 
i=i fc=i 

Here P(fc) is the density of nodes with degree k ^ 
and /(ct, ka) represents the contribution from a single 
fc-degree node. Note that /(ct, ka) actually depends on 
a and ka 



/(ct, fcCT) = /(ct, fccr). 



(40) 



Therefore, we further replace the sum over k in (|39p by an 
integral over k ([3]), and introduce y = ka as the variable 
of integration. 



Tr,(...) = / dslS(L-\sl\)(...). 



(33) 



The ^-function ensures that all spins si have the same 
absolute value L. Substituting (f3Tj) into (f32|) and taking 
the trace (some details of the calculations are given in 
Appendix [Bj one arrives at the free energy per site: 



F{a,T) = -T/N\nZ, 



AIF- 



(34) 



As usual in the mean field approach, the free energy (|34|) 
depends on the macroscopic mean field variable ct. The 
last is to be eliminated by corresponding minimization 
of F(a,T). The expression for the free energy per site 
reads: 



1 N 



(35) 



with 



F(ct,T) =ct 



A-l 



P(y)f(a,y)dy. 



(41) 



The convergence of the integral (|4"Tj) for large y can be de- 
rived from the asymptotic behavior of the function I u {z) 
El: 



(42) 



Namely, substituting (|42[) into (j36|) one finally arrives at 

f(#,y)~V, y->oo. (43) 

The last estimate, together with the power law behavior 
(fT|) proves the convergence of the expression for the free 
energy (J35J) for A > 2. 

The behavior of /(ct, y) for small y — > and a small 
anisotropy parameter u/T <C 1 is characterized by the 
smallest term of the Bessel function expansions 22j in 



f(ff,kiS) = -Tln(2 7 rL)+i-^A; i |H 2 



(36) 



Tin (I (kS) 



+6 



uL 4 
T 



+ 



\t\ 4 



and 



(37) 



In ([36| . I n (z) are modified Bessel functions [22j of the 
first kind: 



In{z) 



1 

2tti 



(38) 



It is instructive to observe that in (|36|) the function / 
depends both on ka and on a (via the second and the 



(z 2 /4)" 



z-fcir^ + g + i)' 



(44) 



where T(p) is Euler gamma function. Now, substituting 
(HI} into (f3l)j) one arrives at: 



1. HJL)\. 2 



f(#,V) = fo + 2 J(T y ; T 



-V 



with 



1_ 49^ (^ 
64^ 8 T T 3 
1 JlLf aj + aj 
384 — ^ +••• 



/o = -Tin (2ttL) + -uL 4 



(45) 



(46) 



We will perform the integration in (|41[) using the expan- 
sion of /(ct, y) (US]) and its asymptotics (14"3")) at y — > oo. 
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Those terms of the expansion (|45[) that are well be- 
haved with respect to the integration in (j4"Tj) may be eas- 
ily integrated. These are the terms, in which y M appears 
with fj, < A — 1. The integration of the remainder of the 
series (|45]l (let us denote it as g(a, y)) needs some special 
care. Using the asymptotic behavior of g{<J,y) at small 
and large values of y, the integration is to be performed 
in the same way as for g(x, y) (see Section Til A[) . to obtain 
the free energy as described above. 

To complete the calculations we now pass from the 
average spin a per nearest neighbor to the mean spin (s) 
of a node. Solving the self-consistency equation (f2"6"|) for 
(s) one finds in a linear approximation in a and u 



(47) 



Substituting (|47p into ([45|) one finally obtains the free 
energy density as 



f((^,k(s)) = f + 



2T 2 T 

T , 49 uL A . , a , . A , 

- "s"~ (48) 
+^«».> 4 +<«m... 

Note, that taking into account higher order corrections 
in (|4"T|) does not change the free energy at critical point. 

The expression (|48|) serves as an example for a mi- 
croscopic interpretation of the phenomenological Landau 
free energy $(af, T) ([3]), ([5]). Indeed, the two compo- 
nent order parameter x in ([5]) may be interpreted as the 
two-component mean spin (magnetization) per site (s) in 
(|48|) . The remaining phenomenological Landau param- 
eters may be found by direct comparison of the expres- 
sion and (|4"5)) . In this way, the value of /o in (T5J) has 
a microscopic representation in terms of (|46p . while the 
coefficients a,- read: 



a = 0, ai = 



2T 2 



J(k) 2 L 4: 



0.2 



(k) 2 L 2 ' 



(49) 



Recall, that in the frames of the Landau approach the 
assumption ao = leads to the dependence of T c on 
(k 2 )/(k). Now we find the expression for T c in the mi- 
croscopic model as: 



Moreover, the microscopic approach allows us to esti- 
mate the temperature dependence of the proportionality 
coefficients in the free energy expressions. The latter is 
of primary importance for the case 2 < A < 3, when the 
critical temperature diverges. Then 

C,C'~T 2 , D, D ~ T, E, E ~ T°. (53) 

In the following, we pass to a more detailed analysis of 
the Landau free energy ([3]). 



III. PHASE DIAGRAMS 

Having determined the behavior of the free energy ([3]) , 
let us investigate the stable states of the system. The 
latter may be found from the minimization of the free 
energy. The condition of stationarity requires the first 
derivatives of the free energy to vanish 



d$(x,T) =Q d$(x,T) =Q 



dx i 



dxo 



(54) 



The stationarity point is a minimum if both eigenvalues 
of the matrix of second derivatives 



_ a 2 $(f,r) 

OX^OXv 



H,v= 1,2 



(55) 



are positive. This condition may also be written as 
ReO^) > 0, det(w jUV ) > 0, fj,,v = l, 2. 



(56) 



From a physical point of view, the minimum of the free 
energy requires positive isothermal susceptibilities. 

In the following we consider the stable states of the 
system with coupled order parameters for the relevant 
ranges of the exponent A, discussed for the generalized 
Landau free energy. 



A. Case A > 5 

For A > 5 the system is described by the Landau free 
energy (TT5)) . whereas the type of the ordering below T c 
depends of the interplay between the fourth-order cou- 
plings. If > and £>W > 0, the system is character- 
ized by order parameter components 



T r . = 



JL 2 (k 2 ) 
2 (k) 



(50) 



The values for the other coefficients in the Landau ex- 
pansion are as follows 



n = xl w (T c -Tf, x 2 = 0. (57) 

If < and 4&( A ) + > 0, both order parameters 
have the same value 



bi = Ci = 0, i = . . . 3, 



4(fc) 4 L' 



■(1 



49 uL 4 



c 4 



24(ky- 



(51) 



(52) 



X\ = x 2 



2 a 



46(A) _|_ C (A) 



(T c - T) 







(58) 



with [3 — 1/2. Here and below we do not write explicitly 
one more solution x\ = 0, x 2 7^ which is symmetric 
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4h) 





^8 




4b (A ^0 





J*) 



FIG. 1: The phase diagram for the coupled two component 
order parameter model (2) on a complex scale- free network. 
The picture shows what type of ordering is realized in the 
different phases depending on the free energy parameters 6^' 
and c (x \ The blank part of the phase diagram corresponds 
to absence of a stable phase. An ordered phase exists only 
if & (A) > 0. The sign of the coefficient c' A ' separates two 
phases. Namely, positive values of c' A ' > correspond to 
phases with only one non-zero order parameter component 
([1, 0] or [0, 1]); negative values c' A ' < that satisfy condition 
(4&' A ' + c' A ' > 0) corresponds to the ordered phase [1, 1], 
where both order parameters to have the same non-zero value. 



to (|57[) and which is stable under the same conditions 
as the solution x\ ^ 0,x 2 = 0. The resulting phase 
diagram is shown in Fig. [TJ The blank parts of the 
phase diagram correspond to cases where no stable state 
exists. For these values of a, and the condition 
of stability of the thermodynamic potential cannot be 
satisfied (i.e. the asymptotics T) — > oo for \x\ — > oo 
do not hold), therefore the system is undefined for this 
range of parameters. 



Note that expressions ([59)1 . ([60)) represent approximate 
solutions of Eq. (f5l)l . which is transcendental when con- 
sidering the free energy (fT8| . Both phases are character- 
ized by the same value of the critical exponent (3 = 1/2. 
For other values of b^' and the free energy (fl"8f does 
not lead to any equilibrium stable state, as in the case 
A > 5. These results are also depicted in the phase dia- 
gram in Fig[TJ 



C. Case 3 < A < 5 

For degree distributions governed by an exponent in 
the range 3 < A < 5, the system is described by the free 
energy (f2"0"|) . Below T c stable states exist only if b^ > 0. 
Namely, there are two stable phases, one with 



Xi 



. (A — 1)&( A ) - 
and a second one with 



x 2 = 0, (61) 



Xl = x 2 



16a 



(To 



Tf 
(62) 

where (5 = . The regions where these states are real- 
ized are shown in the phase diagram FigQ] If > the 
system is in the stable state (|6"Tj) . Otherwise for negative 
c( A ' < and 4&( A ) + > the system is described by 
the order parameter components (|62p . As observed ear- 
lier for larger values of A, the stability conditions (|54p . 
(|56|) cannot be satisfied for other values of b^ and c' A ) 
parameters in the free energy (|20|) . 

We conclude that if (fc 4 ) diverges but (k 2 ) is finite 
(3 < A < 5), the critical behavior differs from the clas- 
sical mean-field behavior. Furthermore, if A = 5, loga- 
rithmic corrections appear, and if 3 < A < 5, the critical 
exponents are functions of A. Note, that for all values 
of A > 3 considered above there exists a finite critical 
temperature. This will not be the case for the values of 
A considered below. 



B. Case A 



D. 



Case 2 < A < 3 



If the exponent A is at its marginal value A = 5, the 
free energy is described by (TT5)) . For temperatures below 
T c stable states exist only if 6 (A) > 0. For c (A) > 0, the 
system is described near the critical point T — > T c by an 
ordered phase with 



2a 



(T c - TY 



bW [IniTc-T)- 1 } 1 / 2 ' 



x 2 = 0. 



If -4&( A ) < c( A) < 0, the ordered phase at T 
characterized by the order parameters 



xi = x 2 



a (T c - T f 

46(A) +C (A) [ln^-T)" 1 ] 1 ^ 



(59) 



T r iS 



(60) 



When the exponent A is in the range 2 < A < 3 and the 
second moment (k 2 ) of the node degree distribution ^ 
diverges, one may infer from ([7]) that the order-disorder 
phase transition does not occur at any finite temperature. 
Taking into account that at T = the system is ordered, 
the system keeps order at any finite temperature, as has 
been confirmed for the Ising model on the infinite size 
Barabasi- Albert scale- free network 0, 0] • 

In the case A = 3, the free energy is given by Eq. (f2l"j) . 
If the parameter D is positive as follows from Eqs. |T]), 
(|23|) . (|49|). the system is always ordered. The type of 
order found depends on the parameter E. If E is positive 
only one order parameter has nonzero value 



Xl 



2C + D 
2D 



x 2 = 0. 



(63) 



8 



For negative values of E both order parameters are 
nonzero and have equal value 



X\ = x 2 



1 _ 4C+2D+E 
f> 4D 

V2 B 



(64) 



The high-temperature dependence of both order param- 
eters in view of (l53l) follows 



Xl,X2 



-vfT 



(65) 



with some constant r\ determined by the coefficients of 
the high temperature behavior of C and D in Eq. (|53|) . 

In the case 2 < A < 3 the system is described by the 
free energy (|22[) . Assuming the anisotropy parameter 
to be small (E <C D ) one finds stable states of the 
system. If C' > (corresponding to iQ, (0SJ) 
and D < 0, the system is always ordered and the type 
of order depends on the sign of the anisotropy parameter 
E . Namely, if E > 0, the ordered phase is characterized 
by 



Xi 



A - 1 



D' 



x 2 = 0. 



If E < 0, both order parameters are nonzero with 

d 



Xl = x 2 



i 

2 2 \ x=3 



A - 1 



AD' + E' 



(66) 



(67) 



Taking into account the high temperature dependence 
of C and D (f53")> . the temperature dependencies of the 
non-zero order parameters for T — ► oo can be found as: 



(68) 



This corresponds to the obtained for the scalar theory 
results [1, 0. As one may expect, for all 2 < A < 3 both 
xi and x-i vanish only at infinitely large temperature. 



IV. REACTION OF THE SYSTEM TO AN 
EXTERNAL ACTION 

A. Isothermal susceptibilities 

In the case of two order parameters, the behavior of the 
system in an external field is described by two quantities. 
The longitudinal susceptibility \\\ describes the reaction 
of the system to the external field applied along the order 
parameter direction. In turn, x± describes the reaction 
to a transverse external field. 

In the disordered state and in absence of an exter- 
nal field the system is isotropic and therefore there is no 
difference between x\\ an d X-L- I n the general case the 
susceptibility matrix x M „ = (dx M / dh v )\T (see e.g. Q) 



(69) 



Here 5^ u is the Kronecker 



depends on both x\\ an d X-L which may be found as the 
inverse eigenvalues of the matrix of second order deriva- 
tives of the free energy ([55 
symbol. 

Thus, above the critical temperature T > T c both sus- 
ceptibilities have the same dependence, for all values of 
A > 3 



1 



A' 1 1 



X-i 



(T-T c 



(70) 



with the standard mean field critical exponent 7 = 1. As 
mentioned above, there is no disordered state in a scale 
free network of infinite size with 2 < A < 3. 

For all A > 3 and below the critical temperature both 
X|| and x± follow a power law with the mean field criti- 
cal exponent 7 = 1, as also found for the case of a scalar 
order parameter. Furthermore, the absolute value of the 
longitudinal susceptibility x\\ coincides with the suscep- 
tibility x found for the scalar case [15j. Our results are 



All 



A > 5 

3 < A < 5. 



(71) 



The absolute value of the transverse susceptibility x± 
depends on both A and the type of order. So, for A > 5 



2?/ 



X-L 



4b<*> 



,(A) 



2ac 



(T c -T)-t, 



for x 
for x 



[1.0] 
[1-1]. 



(72) 



As one may see from (f72"]) . when the coefficient = 0, 
and thus the system described by the free energy © 
becomes isotropic, the transverse susceptibility diverges 
Xi_ — * 00 for any T < T c . This behavior of X-L is quite 
physical and is a consequence of the free energy symme- 
try: an infinitely small external field applied in a direc- 
tion perpendicular to the order parameter, immediately 
changes the order parameter orientation. This is the 
Goldstone phenomenon, corresponding to the existence 
of a soft excitation mode in the ordered phase [23|. 
For 3 < A < 5, the transverse susceptibility is given by 



f (A-1)6< X > rrp _ T N- 7 
Xl - - ) (A-l)(4b< A >+c< A >) 



for x= [1,0] 
^-x (T c -T)-t, for 2 =[1,1] 



(73) 

As discussed above for the case A > 5, again the trans- 
verse susceptibility diverges for a vanishing parameter 
=0. 

For A = 3 the longitudinal susceptibility reads 



1 T~ 



(74) 



The sign of the transverse susceptibilities depends on the 
phase: 



X-L = 



1/2E, 
-1/2E, 



[1,0] 
[1,1]. 



(75) 
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For 2 < A < 3 the behavior is similar. The longitudinal 
susceptibility follows 



1 



XII 



2(3 -\)C 



(76) 



while we have different transverse susceptibility in the 
two stable phases 



The amplitude ratios for the different ranges of A and 
phases are summarized in Table HI 

Summarizing, we note that for all the range of A > 3 
(where the critical temperature T c exists), the behavior of 
the system with respect to an external field is governed by 
a mean- field critical exponent 7 = 1, but the amplitude 
ratios have nontrivial forms. 



X± 



A-l D 
4 q'E' ; 
A-l 4_D +E 



16 



C E 



[1,0] 

[1,1]. 



(77) 



As we learn from the above equations ([70)) - (|75|) . the 
singularity at the critical point is governed by the mean- 
field value of the critical exponent 7 = 1. This repro- 
duces the behavior observed within the Landau theory 
for systems with a scalar order parameter on scale-free 
networks [l5[ . In this respect passing to a system with a 
more complicated symmetry does not appear to modify 
the universal critical exponents. Note however, the sig- 
nificant change in other universal quantities, namely the 
susceptibility amplitude ratios. Defining the amplitudes 
for the susceptibilities by 



Xi 



r-,i(r c 



T c y 

-T)~ 



T > T c 
T<T C . 



(78) 



let us compare the amplitude ratios r + /T_ for longitu- 
dinal and transverse susceptibilities for different phases. 
Recall that for a scalar order parameter Landau theory 
gives r + /r_ = 2 [2(J. Correspondingly, for the free en- 
ergy @ one finds for the longitudinal susceptibility 



(r+/r_)n =2 



(79) 



while the amplitude ratio for the transverse susceptibility 
depends on the type of the ordered phase: 



(r+/r_)_L 



c/2b, £=[1,0] 
-2c/(46 + c), x= [1,1], 



(80) 



where the notations [1,0] and [1,1] indicate the corre- 
sponding phases. As one can see, the amplitude ratios 
([80)) depend on the couplings b, c. For A > 5 the free 
energy (fT5|) is equivalent to that of Eq. @ however with 
coefficients b and c given by Eq. (|17p . Thus the ratio 
r + /r_ attains the same values as for the systems with 
free energy ([2]). For A < 5 the ratio r + /r_ is a function 
of A, similar as it holds for the order parameter critical 
exponent /3. So, the amplitude ratio for the longitudinal 
susceptibility for all 3 < A < 5 reads 



(r+/r_)u = (A-3). 



(81) 



For the transverse susceptibilities the ratio depends on 
the phase, and respectively, on the values of the coeffi- 
cients of the free energy function. Namely, the suscepti- 
bility ratios are 



(r+/r_)j 



2c< A > 



8c' A > 

(A-l)(46( A )+c( A >) ' 



[1,0] 

[1,1] 



(82) 



B. Heat capacity 

The heat capacity describes the behavior of the system 
with respect to a change in temperature 



Ch 



t(-) ■ 

\dTJh 



(83) 



In the frames of the Landau theory, the coefficient of 
I a; 1 2 in the free energy changes its sign at the critical 
temperature. The other coefficients are assumed to be 
temperature independent. Note, that for a phase transi- 
tion on a scale-free network, this assumption holds also 
for A > 3, whereas for 2 < A < 3 the temperature depen- 
dencies of the coefficients are described by (|53[) . Then 
one may find the entropy of the system as 



S 



,9£ 
"<9T 



which for A > 3 reduces to the simple expression 



(84) 



(85) 



where |x| is a function of temperature and external field. 
Substituting stable solutions that follow from (I54|) into 
([84| , one finds the entropy S at fixed external field h = 
for each phase. Respectively, the heat capacity may be 
found by taking the derivative of the entropy with respect 
to the temperature in ([83]) . 

It is known that for a second-order phase transition 
in simple magnets the Landau theory predicts a step in 
the heat capacity at the critical temperature T c . The 
behavior of the heat capacity for a system on a scale- 
free network is richer. In the standard mean-field region 
A > 5 and below T c the heat capacity decreases linearly 
with the decrease of temperature. At the critical point 
the step in the heat capacity is 



Sc h 



2bW 



(86) 



Taking into account the microscopic relations (|4"9"]l . (|50[) . 
(l52)l . one obtains the step in the heat capacity as follows: 



5c h 



2\2 



(87) 



which vanishes for A — > 5. For 3 < A < 5 there is no step 
of Ch at the critical point. Namely, for A = 5 we find the 
following expressions for the heat capacity at T < T c in 
the phases [1,0] and [1,1], correspondingly: 
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A 


(r+/r_)n 


(r+/r-) x[1 , 0] 


(r+/r_) ±[M1 


A > 5 


2 


C ( A '/2& (A) 


-2 C < A V(4fc (A) +c (A) ) 


3 < A < 5 


A-3 


2c( A V(A - l)b (A) 


-8c< A V(A- 1)(4&< A > +c (A) ) 



TABLE I: Amplitude ratios for different ranges of the A exponent. Second column: amplitude ratio for the longitudinal 
susceptibilities; third and fouth columns - for the transverse susceptibilities for two different phases, denoted by [1, 0] and [1, 1] 
respectively. 



Ch 



T 



bW ln(T c -T)- 1 ' 



in phase [1,0], 



Ch 



Aa 1 



T 



4bW + cW ln(T c - T)- 1 ' 
In the case 3 < A < 5 the corresponding formulas read 



in phase [1,1]. 



Ch 



A-3 



4n 



Ch 



A-3 



(A - l)bW 
16a 



2/(A-3) 



(A - 1)(4M A ) +cW) 



T{T c -T)^ 5 -^^ x - 3 \ in phase [1,0], 

T(T c -T)( 5 - A '/( A - 3 \ in phase [1,1]. 
I 



2/(A-3) 



(89) 



(90) 



(91) 



As one can see from Eqs. (188 | (j9Tj) , the heat capacity 
vanishes as T — > T c which differs from the case A > 5, 
where the corresponding value at T c is given by (|8"6"]) . 
Nevertheless, a maximum of Ch is still present for 3 < 
A < 5. Only now, it is shifted from T c to the temperature 
region T < T c . The low temperature behavior of the heat 
capacity at 3 < A < 5 resembles that for A > 5: Ch ~ T. 
The heat capacity vanishes both at T = and T = T c 
and possesses maximum at an intermediate temperature 
< Tq < T c . For A — 5 this temperature coincides with 
T r whereas for lower values of A we find: 



A — 3 

Th = — — T r , 3 < A < 5. 



(92) 



Taking into account the explicit calculations of Section 

MM 

JL 2 

T = (A-2) — , 3<A<5. (93) 

In Fig. [2] we show the typical behavior of Ch for differ- 
ent values of 3 < A < 5. There, we represent Eqs. 
(IfjTjl in the form 



Ch 



T ( 



T\ (5-A)/(A-3) 

t'J 



(94) 



and plot Ch/cQ as a function of a scaled variable T/T c . 

As A approaches from above 3, the critical temperature 
increases and becomes infinite for 2 < A < 3: the system 



0,6 



O 

o 



O 0,4 



0,0 





1 1 


i 1 










/\ k=4.8-' 














X=3.5 


■ ' X=4 \ 




^1=3.1 
— -_i i_ 







T/T r 



FIG. 2: Typical behavior of the heat capacity for different 
values of A in the range of 3 < A < 5. A dotted curve shows 
position of a maximum at temperature To, see Eq. (|92|l . 



is always ordered and the type of ordered phase is gov- 
erned by signs of the coefficients E, E in the Landau free 
energies (|2"TT) . (|2"2"j) . For both ordered phases we obtain 
that the high-temperature behavior of Ch is described by 



(95) 
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where £ depends on the coefficients defined in Eq. (|53|) . 

As we have observed for the order parameter and the 
susceptibility, the character of the temperature depen- 
dence of the heat capacity for the system of two coupled 
order parameters reproduces the one obtained for a single 
scalar order parameter [tj [l5| . Note however the differ- 
ent amplitudes for this behavior resulting from Eqs. (|50")) . 

USD. 



V. CONCLUSIONS 

Models that display phase transitions with two cou- 
pled order parameters serve as archetypes to describe the 
phase behavior in systems with several possible types of 
ordering [l|| [l7|. For example, a system may display 
both ferromagnetic and antiferromagnetic order with a 
coupling between the two. Others may show magnetic 
and superconducting, ferroelectric and ferromagnetic, or 
structural and magnetic order. These models are known 
for their rich phase diagrams and non-trivial critical be- 
havior. Inspired by these observations, this paper inves- 
tigates the phase transitions of a corresponding model 
defined on a scale-free network. Besides the academic in- 
terest, this problem may have implications for models of 
opinion formation on social networks when opinions on 
different issues may be coupled, as e.g. t he p references 
for both a candidate and a political party 18| . 

Our analysis was based on thermodynamic arguments 
in the spirit of Landau theory, as suited for the descrip- 
tion of phase transitions on scale free networks [l5j]. To 
add a microscopic background to the phenomenological 
approach we have also studied a particular spin Hamilto- 
nian that leads to coupled scalar order behavior using the 
mean field approximation. Our results show that for the 
scale free networks with a degree distribution governed 
by an exponent A > 2 the system is characterized by ei- 
ther of two types of ordering. Either one of the two order 
parameters is zero (the [1,0] or the [0, 1] phase) or both 
are non-zero but have the same value (the [1, 1] phase). 
Along with the critical behavior of the scalar order pa- 
rameter systems on scale-free networks, the order of the 
phase transition in the coupled scalar order parameter 
system depends on the strength of the node-degree distri- 
bution decay. For rapidly decaying distributions (A > 5) 



the second-order phase transition is similar to that de- 
scribed by usual Landau theory. Nevertheless the new 
features appear as A decreases: whereas the magnetic 
susceptibility (and higher than the second derivatives of 
the free energy with respect to the magnetic field) remain 
divergent at T c for all 3 < A < 5, the order of the lowest 
divergent temperature derivative of the free energy de- 
pends on A [241 ) . Namely, as seen from Eqs. (fil)0"|) - (p?Tj) . it 
is the third order derivative for 4 < A < 5, the fourth or- 
der for 3 1 < A < 5, and so on until it is only the infinite 
order derivative that diverges for A = 3: the order of the 
phase transition becomes infinite 0, [2l[ . 

The critical behavior of the model considered gives rise 
to non-trivial critical exponents, amplitude ratios and 
susceptibilities. While the critical exponents do not dif- 
fer from those of a model with a single order parameter on 
a scale free network [15J there are notable differences for 
the amplitude ratios and susceptibilities. Another pecu- 
liarity of the model is that the transverse susceptibility is 
divergent at all T < T c ,when O(n) symmetry is present. 
Such behavior is related to the appearance of Goldstone 
modes. It is worth to mention a peculiarity in the be- 
havior of the specific heat. Whereas for A > 5 it has a 
step at the critical temperature T c , this step disappears 
for A < 5. The heat capacity vanishes both at T = 
and T = T c and possesses maximum at an intermediate 
temperature < Tq < T c . 

The phenomena observed serve as evidence of a rich 
critical behavior caused by scale-free properties of the 
underlying network structure. An attractive feature for 
the theoretical analysis of this behavior is of course that 
non-trivial effects are found already in very simple ap- 
proximations. Natural continuations of our study will 
include extensions beyond the mean field approach tak- 
ing into account order parameter fluctuations, further, 
studies of dynamic processes and in particular the criti- 
cal dynamics resulting at or near the critical point. Such 
studies need to be based on more detailed information 
about the structure of the network than the degree dis- 
tribution, such as provided by the adjacency matrix or 
the network Laplacian (Kirchhoff matrix) e.g. in terms 
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of their respective eigenvalue spectra 
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h(s,x, A) 



Sl dy 



(A.4) 



From the asymptotic behavior (1A.2|) one may infer that 
Qi is convergent. Assuming that near the critical point 
the absolute value of the order parameter |a;| <C 1 is small, 
we replace the function g(x,y) in Q 2 by its expansion 
(jA.ljl at small values of y and perform the integration. 
Then one obtains 



dy 



X-\ ] ,To \ \X 



y 

where 



—g(x, y) = Q(x, A) - (& 4 + c 4 1 2 ) 



5-A 



5- A 
(A.5) 



xi + xi, £ 5 ~ A 



Q(x,X) = Qi(e,x,X) + (b 4 + Ci 1 2 ) 



5-A 



(A.6) 



Naturally, Q(x, A) does not depend on e (as Jj2" ^g(x, y) 
does not depend on e), while the dependence of Q(x, A) 
on x is reasonably (see the asymptotics (|A. 1|) and 
Eq. (|A.6|) ) to be assumed as follows: 



4 , 4 

/ -* w X-t ~r Xn 

Q{x, A) = vi + v 2 4 



(A.7) 



where V\ and t>2 are some coefficients, in general depen- 
dent on A and the temperature. 

Substituting these results into (|13|) . one obtains: 



dkP{k)g{x,k\x\) = AQ(x,X)\x 



A 



A-l 



+ ^ + c *%^' 4 - < A - 



In the region of A considered (4 < A < 5) near the crit- 
ical point the leading term includes a factor |cc| A ~ 1 and 
correspondingly Q(x, A) is part of the relevant terms of 
the free energy. 



APPENDIX A 



APPENDIX B 



In order to perform the integration in (fTI]) . assume 
4 < A < 5. For other values of the exponent A the 
calculations can be performed in a similar way. From 
Eqs. flU)) and (fl"2"|) we derive the following asymptotics of 

g(x,y)- 



Here we calculate the partition function (|32|) with 
H % MF described by (f3"Tj) . To calculate 



N 



N 



Zmf = J] Z mf = II Tr * e~ H ' MF/T (B.l) 



g(x,y) 

g{x,y) 



(b4 



■ c 4 - 

v - 



\xr 

oo . 



y^0 



To analyze Eq. (|T4|) let us define 

Qi(s,x,X) = / —g(x,y), 
J e y 



(A.l) 
(A.2) 

(A.3) 



we use the following property of the (5-function 

( 5(L-| S 1|) = 2L,5( J L 2 -| S 1| 2 ) 
and use its Fourier presentation 



1 / >io< 
S(x) = — \ 



duoe 1 



(B.2) 



(B.3) 



13 



Then Z l MF reads 



to obtain 



LZi 



oc />100 



Z MF = —r / dsl / dv e 



„v L 2 



(B.4) 



oo J — 100 



with 



Now, let us use the representation 

exp{-M 4,i - v o s l,i + k°uS v ,i} 

= expl-^^-jexpl-wos^i + ji^s^i} (B.6) 

where exp(— |f ^4) is interpreted as 



exp{- 



} = 1 



Uq <9 4 1 /Mo 5 



+ +■•• (B.7) 



j 4 dcr 4j jfdat ' 21" jfdat 

Substituting (|B.6|) into (|B.4|) one obtains 

71-1 V ,7=1 3t da » ■/-«= 

/ioo 2 
dv a e VoL2 Y[ exp{-«os* 4 + 7,1,--,,., } (B.8) 
- io ° ^=1 

To change the order of integration over S Vj i and i>o , we 
multiply the integrand with exp{a(L 2 — |si| 2 )}, which is 
equal to unity due to the constraint. Let us choose a to 
be sufficiently large to satisfy (vq + o)si li — jiO v s v ,i > 0. 
Then one may use the Poisson integral 



dxe~ ax +bx = J-e b /4a (B.9) 
a 



^MF — 



Lzi 



a+ioc 



dve 



uL- 



n 

i/=l 



e If s^FgjTo-;/^ 



(B.10) 

Assuming the anisotropy parameter u to be small, and 
respectively uq <C 1, we keep only the term linear in uq 
in the (|B.10p expansion. Then (|B.10j) may be written as 



nri 



XUJ- 



^ duJi exp{ ^jtLaiuJi + 1/w,)} 



l-u Q L 4 {6 



6- 



where 



{jiLa) 2 j l La 



2Lv 



(B.ll) 



(B.12) 



The integration path for the variable u>i in the integral 
in (jB.lip ranges from 2La/ jiO — ioo to 2La/ji<j + ioo. 

Using the definition of the modified Bessel function of 
the first kind 



In(z) = 

one may write 



J MF 



1 

2nd 



2-kLz 



e ( 2: /2)(a;+l/c) a; -n-l dW) 



(B.13) 



h{jiLu) - u L 4 {6 



{jiLcr? 



h(jiLa) erf + cr| , 



(B.14) 



Substituting (|B~T4f into (|B~1~1) one finally obtains the par- 
tition function, and respectively one may find the free en- 
ergy per site F(a,T) = —T/NXyiZmf- Again, keeping 
terms linear in u, one obtains the free energy (|35D . 



